A local Hamiltonian for spherically symmetric gravity coupled to a scalar field 
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We present a gauge fixing of gravity coupled to a scalar field in spherical symmetry such that 
the Hamiltonian is an integral over space of a local density. Such a formulation had proved elusive 
over the years. As in any gauge fixing, it works for a restricted set of initial data. We argue that 
the set could be large enough to attempt a quantization the could include the important case of an 
evaporating black hole. 
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Spherically symmetric gravity coupled to a scalar field has been an arena where many seminal ideas of black hole 
physics originated, through classical and semi-classical treatments. The full quantization of the model has resisted 
analysis, in part due to the complexity of the Hamiltonian structure of the system. If one could gain control of this 
model it would be a superb scenario to test key ideas about black hole evaporation. The first attempt to treat the 
problem quantum mechanically was carried out by Berger, Chitre, Nutku and Moncrief [l[ and further developed by 
Unruh [2| • The resulting Hamiltonian was intractable enough that Unruh remarked "I present it here in the hope that 
someone else may be able to do something with it." More recently, Husain and Winkler and Daghigh, Kunstatter 
and Gegenberg Q, using Painleve-Gullstrand coordinates simplified somewhat Unruh's treatment. All these efforts, 
however, failed to produce a local Hamiltonian. We would like to show that using Ashtekar's new variables a gauge 
fixing can be found that yields a Hamiltonian that is the spatial integral of a Hamiltonian density. A similar gauge 
fixing can be carried out in traditional variables [3| ■ It also appears to apply in other 1 + 1 models, like the Callan- 
Giddings-Strominger-Horowitz black holes [5(- We do not have a clear explanation as to why it seems to apply in 
such generality, it appears to be related to the possibility of defining a mass function [H, 0] • 

The subject of spherical symmetry with Ashtekar's new variables has been discussed in many instances. We will 
not carry out a full discussion here. We refer the readers to the literature. This is just a minimal introduction in 
order to make the paper self-consistent. The topology of the spatial manifold will be chosen of the form S = R + x S 2 . 
We will use a radial coordinate x and study the theory in the range [0, oo]. The case in which there is a horizon at 
x = can be treated with suitable boundary conditions. 

The formalism for dealing with spherically symmetric gravity with Ashtekar's new variables was discussed by 
Bojowald and Swiderski[8| and also in our recent paper[9(. It is best to make several changes of variables to simplify 
things and improve asymptotic behaviors. We will not go through all these steps here. It suffices to notice that at 
the end of the process one is left with two pair of canonical variables E v and K v (in our recent paper called A v ), 
and E x and K x that are related to the traditional canonical variables in spherical symmetry 



ds 2 = A 2 dx 2 + R 2 dn 2 
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and Pa = V E x K lp /(2 r y) where 7 is the Barbero-Immirzi parameter and Pa is the momentum canonically conjugate 
to A and we are considering the positive branch of E x . One also has that the conjugate momentum to the variable 



R is given by Pr 



-. The diffeomorphism and Hamiltonian constraint can be written as 1 , 
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Recalling that the total Hamiltonian for the system is given by Ht = J dx{N x C x + NH), one can redefine the 
shift N x cw = A^ ld + 2NK v y/E~ x / (E x )\ and the lapse A^ new = Af old (E x )' / E v , one gets a Hamiltonian constraint that 



1 To derive equation J3} one substitutes A v = 2^K V , 2^K X = A x + r)' in equation (47) of reference [i| and adding the scalar field 
contribution, e.g equation (2) of |10| . We have absorbed a factor of 4-7T in Newton's constant. 
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The quantity in the square bracket above is a total derivative, 
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This remarkable property is the key element in allowing to define a local Hamiltonian. Choosing a gauge in which the 
term involving the derivative does not depend on the gravitational variables, one is left with a Hamiltonian that only 
depends algebraically on the gravitational variables. As we mentioned, it appears that this is typical of all theories 
in 1 + 1 dimensions that involve a mass function. It at least holds spherically symmetric Einstein gravity in the 
traditional and new variables and for the Callan-Giddings-Harvey-Strominger model. 
We will completely gauge fix the theory. The first gauge condition is xi — with 
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E x 
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In order to preserve the constraint in time the Lagrange multiplier N x gets fixed N x = The diffcomorphism 
constraint can be solved, determining the variable K x . The only constraint left is the Hamiltonian, which (omitting 
an overall factor l/(G(E'- p ) 2 ) becomes, 
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Our strategy will be to perform a canonical transformation from the variables <j>, Pq, K v , E v to a new set of variables 
X, Px , /, Pf such that X is essentially what appears in the square bracket differentiated. We will later fix the gauge 
by setting X equal to a given function of t,x. As a consequence Px, the canonical momentum of X, will not appear 
differentiated in the constraint. This means that preserving the gauge fixing condition will lead to an algebraic 
equation that determines the lapse, and therefore to a local true Hamiltonian. 

To construct the canonical transformation, let us start by identifying the variable X, 
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Recalling that the scalar field has dimensions of inverse length in 3 + 1 dimensions, the factor x 2 on the right is 
chosen so X has dimensions of length (or time, since we chose c = K = 1), since it will later play the role of time. 
The factor <fi 2 is chosen so weak fields behave well in the gauge fixing (for instance if <f> — one has K v = and 
E^ = xrf\-2GMj x in the usual Schwarzschild gauge). We added a function of time M(t). Later, if one studies 
the fields asymptotically one finds that M is a constant that corresponds to the ADM mass. At the moment it is just 
a choice in the definition of X . 

To complete the canonical transformation we then seek a generating function, we choose it to be of type I 
Fi(<j>, K v , X, /), so one has that (recalling that {K v (x), E v (y)} = G5(x — y)), 
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We start from the first equation and note that we can use ((U) to write in terms of the quantities that the 
generating function depends on, 

E* = ~, (13) 
where we chose the minus sign of the square root so the Hamiltonian is positive definite and for brevity we write, 



Y = dK% + i- ^M. - x X(f>2. (14) 

So we can now proceed to integrate © and choosing the integration constant to give the simplest form to the 
generating function yields, F\ = — |j log (K v + Y) + cf)f . 

With the generating function and (|10j|12l) we find the explicit form of the new variables in terms of the old ones, 

x 2 X<b , 

/, (15) 
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x 2 <\> 2 
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The last equation will become the Hamiltonian constraint when we rewrite the right hand side entirely in terms of 
the new variables. Rewriting <fi is immediate. To obtain K v we solve (JTJ) rewritten in terms of the new variables, i.e., 



(x 2 <j) 2 X)' 2Gx 2 K s „4>'P 
Y 2 + Y 

and P^, and <p given by (|15j) and (|16|) respectively, so we have 



H = + ^ ^ + GP 2 + Gx 4 , (18) 



xU 2 - x + 2MG + XP 2 x 2 



K « ^ 2xU ' ~ ' (19) 



with 

U = 

and 



-. ^ ((P') 2 x 4 GZ + 2x 5 P'XP f + 2 (XP 2 x 2 )' x 3 + 2Xx 3 P f f + f 2 GZ + 2xV?) 2 , (20) 

(/ + x 2 P' f ) 4x~G V ' 



V = x 



(P' f ) 2 xGZ+(XP 2 x 2 )') (f 2 GZ + 2x 3 XP f (x 2 P' f + f) + (XP 2 x 2 )' x 3 ) (21) 

Z = x 2 XP 2 + 2GM — x. (22) 
We now consider (|17|) written entirely in terms of the new variables, 

x 3 P 2 

Px + 2GY{J,P f ) (K v (f,P f ) + Y(f,P f )) = °- (23) 

This expression is the Hamiltonian constraint that is now easy to deparameterize. The total Hamiltonian is given by, 
-ffibtai = J dxN (Px + "Htiuc), where we recognize the true Hamiltonian density, 

x 3 P 2 

%TtnC = 2GY(f,P f )(KM,Pf) + Y(f,P f ))- (24) 

To prove that indeed this expression is the true Hamiltonian density, we proceed to completely fix the gauge. We 
choose X2 = —X + g(x) +t = 0. The preservation in time of this condition, + {^-^Totai} = implies that the 
lapse N = 1. The system is now totally described in terms of the matter field variables /, Pf, since X is fixed by the 
gauge fixing and Px is given by minus the true Hamiltonian. If we now consider the time evolution of the remaining 
variables, 

f = {f,H Total } = {f,H True }, (25) 
P f = {P f , H Total } = {P f , H Tluc } , (26) 
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showing that the true Hamiltonian indeed generates the evolution. 

The expression for K v (fHj| contains a series of square roots. This reflects the fact that the construction will not 
work for generic initial data, as one expects in gauge fixed treatments. In order to analyze under which conditions 
the construction works, we study the situation of weak fields, so we will assume / = O(e) and Pf = 0{e) with e«l 
and we will keep only leading terms in e in all equations. We will also assume that M\fG ^> 1 (we are using units 
where h = c = 1). In order to simplify expressions we will also assume g(x) = cx with c a positive constant. The 
expression for U becomes, 

u = ttW vI iy^' + cxfPf + xZp, f cPf + [( a;3cp /)' { 2cxfp i + + 2x3p f cp i)\ 5 ) 2 • ( 27 ) 

Sufficient conditions for the existence of the square roots are, 

(x 3 cPf)' = w{x), (28) 
(2cxfP f + 2x 3 P' f cP f ) = v(x), (29) 

with w(x) and v(x) positive functions. Solving the differential equations we get, 

\ ix f x w(x')dx' 

Pf = )LJ ~^ , (30) 

3x J x w{x')dx' — x 2 w{x) v(x) 
1 = 2x*P f (x) { ' 

So we see that indeed one can specify initial data in the gauge we chose. 

We have therefore presented for the first time a local Hamiltonian for a scalar field coupled to gravity in spherical 
symmetry, a problem that was unclear had a solution. The technique appears applicable in other 1 + 1 dimensional 
situations where there exists a mass function. The result has a counterpart in path integral treatments, where authors 
were able to integrate out the gravitational variables lll f . This includes the Callan-Giddings-Horowitz-Strominger 
model, which has received renewed attention recently [12| and is one of the best understood models of black hole 
evaporation. In further work we will discuss the boundary treatment in these coordinates and will show the evolution 
of collapsing scalar field pulses numerically. The resulting unconstrained system can be useful for quantization in 
situations involving gravitational collapse and black hole evaporation. 
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